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STOCHASTIC  ERROR  MODELING  OF  BEAMFORMER  OUTPUT 
FOR  ARRAYS  WITH  DIRECTIVE  ELEMENTS 


1.  INTRODUCTION 


The  fundamental  design  goal  of  a  passive  array  is  to  detect  weak  signals  over  a  broad  band  of 
frequencies  in  the  presence  of  noise  that  is  due  to  structures,  water  flow,  the  surrounding  ambient 
environment,  and  extraneous  targets.  The  array  design  is  accomplished  by  matching  some  spatial 
filtering  and  signal  filtering  (commonly  known  as  beamforming)  to  the  expected  noise  and  signal 
characteristics.  The  broadband  beamforming  is  achieved  through  careful  application  of  a  time 
delay  to  each  element  of  the  array  (or  phase  delay  in  the  frequency  domain),  which  creates  a 
passive  spatial/temporal  filter  for  the  expected  signal  and  noise  characteristics.  The  elements  in 
the  final  array  are  expected  to  behave  as  modeled  mathematically  (i.e.,  the  performance  of  the 
array  elements  will  match  the  mathematical  model  used  for  the  beamformer). 

When  arrays  are  manufactured  as  integral  components  of  structures  (as  in  conformal/hull 
arrays),  isolating  the  array  components  from  inherent  manufacturing  uncertainties  becomes  a 
problem  that  can  be  addressed  through  array  tolerance  analysis.  This  problem  has  historically 
been  treated  both  in  textbooks^’^  (general  arrays)  and  in  technical  reports^’^  (sonar  arrays).  In  this 
report,  the  array  tolerance  analysis  is  developed  with  two  new  approaches  to  the  previous 
research. 

The  first  approach  involves  the  addition  of  the  effects  of  directional  sensors  on  the  curved 
surface  of  the  array,  which  is  especially  desirable  for  a  conformal  application  because  the  sensors 
can  then  reject  much  of  the  structural  noise  through  their  directive  response.  However,  the  added 
directivity  makes  the  phased-element  addition  task  of  the  beamformer  very  sensitive  to  certain 
errors. 

The  second  area  addressed  in  this  report  concerns  the  handling  of  failures.  Failures  are 
considered  in  the  traditional  sense,  as  well  as  for  the  case  of  failure  compensation  through 
companion  element  blanking.  Both  techniques  disregard  the  response  from  certain  companion 
elements  (to  the  failed  element)  in  the  array  to  preserve  specific  array  symmetries  in  the  presence 
of  failures.  Two  basic  types  are  considered:  diagonal  element  failure  compensation  and 
beampattem  symmetry  failure  compensation.  Diagonal  element  failure  compensation  is  an 
approach  to  array  element  failures  that  effectively  zeros-out  elements  that  are  geometrically 
diagonal  (relative  to  the  acoustic  center)  to  the  failed  element,  thus  maintaining  the  array’s 
acoustic  center.  Beampattem  symmetry  failure  compensation  effectively  zeros-out  those  elements 
that  are  the  same  distance  from  the  acoustic  center  as  the  failed  sensor,  thus  maintaining  the  array 
symmetry  in  either  the  horizontal  or  vertical  beampattem.  Other  failure  compensation  techniques 
are  not  treated  here  because  they  are  beyond  the  scope  of  this  effort. 
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For  the  cases  considered,  it  is  assumed  that  deterministic  effects  not  dependent  on  statistical 
parameters  are  included  in  the  governing  deterministic  beamforming  expression,  where  they  are 
evaluated  computationally.  The  required  statistical  array  tolerance  effects  are  shown  in  a 
progressive  stochastic  analysis. 

Sections  2  and  3  of  this  report  present  some  mathematical  methods  for  array  beamforming 
and  statistical  analysis  that  are  then  used  throughout  the  analysis.  The  error/tolerance  analysis 
begins  in  section  4  with  the  case  of  array  failures,  both  compensated  and  uncompensated.  Section 
5  presents  the  next  stochastic  addition,  which  treats  element  phase  variations,  and  section  6  uses 
element  amplitude  variations  to  complete  the  stochastic  analysis.  The  important  equations  from 
the  error  analysis  are  summarized  and  some  computational  examples  are  shown  in  section  7. 
Finally,  in  section  8,  examples  of  physically  relevant  errors  and  their  effects  due  to  phase  or 
amplitude  variations  are  described. 
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2.  ARRAY  BEAMFORMER 


A  general  beamforming  expression  in  the  frequency  domain  for  an  array  of  M  elements  is 
given  as 


^  -  r.  -  1 

;)  =  ^  Wm  dm{Q)  P  ’ 


where  h  =  Qa; /c  and  ks  =  are  the  wavevectors  in  the  incident  and  steer  (or  look) 
directions,  respectively.  In  this  notation,  C  represents  a  unit  vector  in  three-dimensional  space 
(hence,  HCll  =  1)-  completeness,  the  time-domain  beamform  output  associated  with  equation 

(1)  is  given  by 


M 

vo{i)  =  S  dm  (Ci)  Um 

m=l 


^  _  (Ci  ~  C)  • 


(2) 


where  Um{t)  is  the  voltage  output  waveform  of  the  m*'"  hydrophone.  These  two  expressions  are 
related  through  a  standard  Fourier  transform. 

The  beamforming  expression  in  equation  (1)  includes  an  element  directivity  terai,  dmiCi)^ 
which  depends  on  the  incident  wave  direction  Q-  For  an  example  of  this  directivity  tem,  it  is 
assumed  that  the  element  directionality  is  given  by  a  cosine  function  with  perfect  baffling  as  m  a 
notional  velocity  sensor;  hence, 

dmiCi)  =  max{Ci  •  rim,  0}, 

where  Um  is  the  unit  normal  to  the  m*"  array  element.  Because  the  dot  product  of  two  unit  vectors 
is  equal  to  the  cosine  of  the  angle  between  them,  this  expression  represents  a  cosine  directionality. 
Also  seen  in  this  equation  is  the  parametric  dependence  of  element  directivity  on  the  element 
normal  (or  orientation).  Thus,  it  is  shown  that  any  variations  in  element  onentation  will  directly 
affect  element  amplitude  through  the  directivity  expression. 
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3.  STATISTICAL  QUANTITIES  OF  INTEREST 


The  general  beamformer  expression  is  written  as  a  sum  of  complex  variables  v^.  in  the  form 

=  S  ^rn  dm{Q  exp  [i  {h  -  k,)  ■  Xm\  , 

m— 1 

=  (4) 

where  the  unmarked  dependencies  are  all  implied.  In  the  context  of  this  analysis,  the  above 
expression  is  a  sum  of  complex  random  variables  Um  that  represent  the  output  of  a  random 
process.  The  quantity  of  interest  when  processing  this  beamformer  output  is  the  beamformer 
output  power  spectrum 

in«)i"  =  (E (E . 

=  (5) 

m  m’ 

where  (•)*  represents  a  complex  conjugate. 

The  statistics  of  are  important  to  the  analysis  of  tolerance  errors  during  array 

manufacturing,  as  well  as  to  the  effect  of  these  errors  on  array  performance.  In  this  study,  the 
statistics  derived  denote  the  expected  value  of  the  random  output  variable.  Throughout  the  sequel, 
the  notation  E{x]  is  used  to  represent  the  expected  value  and  Var{x}  is  used  to  represent  the 
variance  of  a  complex  random  variable  x.  While  the  variance  is  not  derived  explicitly,  the 
techniques  and  simplifications  developed  are  easily  applied  to  derive  expressions  for  variance.  It 
is  important  to  note  that 

E{f(x)}  #  f{E{x}) ,  («) 

with  the  exception  being  for  the  linear  function  /  (x),  so  that,  in  general, 

£{|VM|"}  JS  .  0) 

Thus,  a  simplistic  analysis  provides  no  useful  results. 

For  the  statistics  on  the  beamformer  output  power  spectrum,  the  results  of  Nuttall^  are  applied 
as  shown  below.  The  expected  value  of  the  power  (magnitude  squared)  of  a  sum  of  complex 
variables  gives 


(8) 


=  E  - 1® K>1')  + 
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which  will  be  used  throughout  the  sequel.  For  notational  convenience,  the  implied  dependence  on 
u  is  omitted,  except  for  cases  where  it  serves  to  clarify  the  narrative.  Nuttall  also  shows  that  the 
variance  of  the  power  of  a  sum  of  complex  variables  can  be  written  as 


Var{\V{oj)\‘^} 


£;{|v|‘}  -  (£{|yn)' , 

S  -  (£{|vp})' , 

klmn 


(9) 


where  the  second  term  in  the  last  expression  is  the  square  of  the  previously  derived  mean.  The 
fourth-order  sum  in  this  equation  is  very  complicated  and,  for  that  reason,  cannot  be  written  out 
here  for  the  general  case.  In  evaluating  the  fourth-order  sum,  Nuttall'*  shows  that  only  five 
statistical  quantities  are  required: 


E{v„},  E{\v„W  E{vl},  and  (10) 

Although,  in  the  analysis  that  follows,  only  the  mean  value  is  determined,  the  interested  reader 
can  use  these  techniques  to  determine  values  of  the  above  quantities.  The  values  can  then  be 
applied  to  the  expressions  found  in  Nuttall'*  and  evaluated  numerically.  The  expressions  that 
combine  these  terms  into  the  fourth-order  sum  in  equation  (9)  are  also  given  in  reference  4. 
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4.  ELEMENT  FAILURES 


The  first  stochastic  variation  under  consideration  is  element  failure.  For  array  failure  analysis, 
a  binomial  random  variable  fm  is  used  to  model  the  failure  state  of  the  element.  The  two 
states  of  this  binomial  variable  are  an  element  failure,  fm  =  nnd  a  working  element,  fm  =  0. 

The  failure  variable  impacts  the  term  in  the  beamformer  sum  of  equation  (1)  by  a 
multiplication  factor  of  1-  fm-  The  probability  density  function  or  probability  distribution  of  the 
binomial  random  variable  fm  is  given  by 

P{fm)  =  QSifm  -  1)  +  (1  -  QWm)  ,  (H) 

where  5(-)  is  the  Dirac  delta  and  Q  is  the  individual  element  failure  rate  (i.e.,  an  element  is  failed 
with  probability  Q).  For  the  distribution  given  above,  the  stochastic  moments  of  the  quantity 
1  fm  are  given  by 

E{{1- /„)'•}  =  I  =  Vfc>0,  (12) 

which  provides  all  the  moments  necessary  for  a  stochastic  analysis  of  the  element  failure  for 
independent  failure  events  within  the  beamformer  expression. 

When  element  failure  compensation  techniques  are  applied,  element  failures  become 
statistically  dependent  events.  From  the  beamformer  output  perspective,  this  means  that  elements 
always  fail  in  pairs  (i.e.,  if  one  element  fails,  the  corresponding  zeroed-out  element  does  not 
contribute  to  the  beamformer  output  and  has  thus  “failed”).  Because  these  pairs  are  related 
through  some  geometric  symmetry  (either  left-right,  up-down,  or  diagonal,  depending  on  the 
failure  compensation  scheme),  there  is  a  strong  geometrical  dependence  on  the  failure,  which 
must  be  accounted  for  in  the  mathematical  model. 

To  maintain  generality,  it  is  assumed  that  the  array  elements  are  numbered  such  that  the 
+  \-  m)*'*  element  is  symmetrically  opposing  the  element,  in  the  manner  corresponding 
to  the  element  failure  compensation  scheme  employed.  Now  the  beamforrmng  expression  (1)  is 
rewritten  as  (for  convenience,  it  is  assumed  that  M  is  even  and  odd  M  is  &  simple  extension) 

M/2 

Vo{uj)=  [Z,n  +  ZA]  ,  (13) 

m=l 

with  Zm  representing  the  nominal  beamformer  contribution  of  the  element  as 

Zm  =  Wm  dm{Q  exp  [i  {k  -  kg)  ■  fm]  (14) 

and  m  =  M  -H  1  -  m.  The  failure  of  independent  elements  is  governed  by  the  failure  variable  fm 
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defined  above,  and  accounting  for  compensation,  the  failure  enters  the  beamforming  expression  as 

MI2 

m=l 

M/2 

“  (1  “  /m)(l  /rh)  \^7n  “I"  ^rh]  5  (15) 

m=l 

where  the  subscript  1  on  Vi{u))  is  used  to  represent  the  case  of  element  failures.  This  expression 
shows  that  if  either  the  or  the  element  fails  (either  /^  =  1  or  =  1),  neither  element 
contributes  to  the  beamformer  output 

In  equation  (15),  the  product  of  failures  is  the  only  stochastic  quantity  that  appears.  This 
product  is  denoted  as 

V 

An  ~  (1  ~  /m)(l  ~  /m))  (16) 

such  that 

M/2 

1^1(0;)  =  ^2  0m{Zm.  +  Zm)  ,  (17) 

m=l 

where  the  bivariate  random  variables  {/0m}m=i  independent  of  one  another.  The  stochastic 
moments  of  this  random  variable  are  given  by 

=  II  &P  (Up  (UdUU  ={1-Q)\  V  (:  >  0 .  (18) 

which  is  used  throughout  the  study.  With  these  moments  defined,  the  expected  value  of  the 
beamformer  output  power  spectrum  is  evaluated  as 

M/2 

=  E  (B{lxmR -  +  IBIV'JP,  (19) 

m=:l 

where  . 

Xm  ~  PmiZm  +  Z^h) 

and  a  dependence  on  u  is  assumed. 

The  first  term  in  equation  (19)  is  given  by 

^^{IXmP}  =  E{|/5m(2„  +  2A)P}, 

=  E{el}-\z^  +  z^W 

(1  Q)  ■  \Zm  "h  Zjfi\  ;  (20) 
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the  second  term  is  given  by 


|■E'{Xm}P  —  +  ^m)}l  j 

=  \E{^m}?-\{Zrn^Z^)W 

—  (1  —  Q)^  •  l^m  +  -2mP  ; 

and  the  third  and  final  term  is  given  by 

\E{Vi}f  =  \E{Y^Pm{Zm  +  Z^)}\\ 

=  \J2E{Pm}{Zm  +  Zrn)\‘', 

—  I  ^(1  ~  Q)^(^Tn  +  ^m)  I  7 

=  (i-cr'lE(^".  +  -^)l". 

=  (i-Or-N".  (22) 

Combining  these  terms  back  into  equation  (19)  yields 

M/2 

E{\Vi\^}  =  [(1  -  QY  -  (1  -  Q)1  E  +  (1  -  <5)"^  l^oP  •  (23) 

m—1 

It  is  noted  that  an  element  failure  rate  Q  of  zero  produces 

E{\Vi\'^}  ^  |Fo|2,  as  (3->0,  (24) 

as  expected. 

The  case  of  uncompensated  failures  is  formulated  in  the  same  way  to  arrive  at 

M/2 

=  [(1  -  Q)  -  (1  -  «)"]  E  +  (1  -  ’ 

m=l 

M 

=  [(1  -  Q)  -  (1  -  <3)^1  E  wl  4.(5)  +  (1  -  <3)"  ,  (25) 

m=l 

which  differs  greatly  from  equation  (23).  Thus,  the  use  of  statistics  without  the  geometrical 
considerations  of  compensation  (such  as  doubling  Q)  is  unwise  and  would  yield  highly 
questionable  results.  The  uncompensated  failure  case  is  obtained  from  the  compensated  failure 
case  by  removing  the  cross-element  terms  in  the  sum  of  \zm  +  z^l"^  and  by  replacing 
(1  _  Q)2  (1  _  Q).  This  result  is  used  throughout  the  remaining  derivations. 
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5.  RANDOM  PHASE  VARIATIONS 


The  second  stochastic  variation  under  consideration  is  random  phase.  Such  variations 
represent  either  a  byproduct  of  the  elements  and  their  associated  electronics  or  a  component  of 
random  positional  variations.  The  fonner  is  modeled  as  a  normally  distributed  (Gaussian)  random 
phase  variable  with  zero  mean  and  a  specified  variance.  The  latter  is  modeled  as  a  phase  variation 
of  the  form  (fcj  -  ks)  •  Axm,  where  is  a  Gaussian  random  positional  variation,  with  random 
amplitude  and  direction.  Because  of  the  potential  for  very  complicated  probability  density 
functions  for  the  phase  variation,  a  result  for  general  phase  variations  is  presented.  The  only 
restrictions  on  these  variations  are  that  they  are  statistically  independent  and  identically 
distributed  over  the  elements  m. 

Beginning  with  a  phase  variation  <t>m  that  enters  the  beamforming  expression  for  the  element 
failures  shown  in  equation  (15)  yields 


M/2 

V'2(W)  =tl^rn  ,  (26) 

m=l 

where  the  subscript  2  denotes  element  failures  plus  phase  variations  and  Prn  is  defined  as  before. 
Letting  the  variable  •jm  be  defined  by 


Tm  —  €Xp(2<^77j)  , 


(27) 


the  expression  is  rewritten  as 


M/2 

V2(w)  =  /5m  {‘1/m^m  "b  'lifi^rh)  •  (28) 

m=l 

The  stochastic  moments  of  the  variable  7^  are  given  by 

£{7™}  =  /  ,  k>0,  (29) 

where  p  {<j)m)  is  the  probability  density  function  of  the  random  phase  variable  4>rn-  To  maintain 
generality,  the  only  assumptions  for  the  specific  form  of  the  random  variable  are  that  4)^  and 
<f)n  are  independent  and  identically  distributed  for  all  m^n,  so  that  ^{7^}  =  £'{7*}  for  all 
A;  >  0  and  all  m,  n. 

As  for  the  case  of  element  failure,  the  expected  value  of  the  beamformer  output  power 
spectrum  is  evaluated  as 


M/2 


m=l 


(30) 
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where 


and  a  dependence  on  to  is  assumed. 

To  evaluate  the  terms  in  equation  (30),  two  properties  of  7^  are  derived.  The  first  is  given  by 

bm?  =  7m7m  =  exp(i<^;n)  exp{-i(j)m)  =  1 ,  (31) 

which  is  independent  of  the  random  variable  (pm-  The  second  property  is 

Ebmi:}  =  -E{7.n}S{7;}  =  Eb,.}  (£{7„})-  =  ,  (32) 

which  follows  from  7„  and  7^  being  independent  and  identically  distributed.  The  two  properties, 
shown  in  equations  (31)  and  (32),  are  used  to  simplify  the  expressions  in  equation  (30). 

The  first  term  in  expression  (30)  is  given  by 

-  E{\^m{lmZm+lmZm)?}  , 

~  }■ -^{177711  \Zm\  |7m|  \Zm\  ’ym'Ym^TnZT^  ~i- 'Jm'yrhZmZm}' y 

=  (1  -  Q)2  •  {^Zm?  +  |2mP  +  \E{'^m}?{z*^Zm  +  2m2*))  i  (33) 

the  second  term  is  given  by 

\E{'ipm}?  =  \E{^m{lmZm-\-lmZm)}?  y 

=  \E{^m]?  ■  \E{'imZm  +  ImZ-fh)}?  y 

=  (1  -  Qf  ■  \E{^m}{Zm  +  2m) P  , 

=  {^~QT '\E{'Ym}^  '\Zm-^  Zm\^  \  (34) 

and  the  final  term  is  given  by 

\E{V2}?  =  \E{Y,Pm{lmZm+l^Zm)}\\ 

=  I  I]  E{^m}E{''imZm  +  ImZm}?  y 

=  •{E{lm}Zm  +  E{-f^}Zm)\^  , 

=  \{l-QfE{^m}Y.i^m  +  Zm)\\ 

=  (l-Q)^.l-E{7m}p-ll^oP,  (35) 
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which  makes  use  of  the  independence  of  on  m.  Combining  these  expressions  back  into 

equation  (30)  and  then  performing  some  algebra  yields 


771=1 

M/2 

+  [G2  (1  -  Q)2  -  (1  -  Q)4]  E  (2m4  + 

m=l 


(36) 


+  GH1-  QY  |l/oP, 


where 

G  =  |E{7m}l  =  |B{e‘*"}| 


is  written  to  illustrate  the  independence  on  m. 

Three  notable  limits  are  found  in  equation  (36).  First,  when  there  are  no  phase  variations,  the 
expression  becomes 

M/2 

E{\V2\^}  [{I-QY  -{l-QY]i,\Zm  +  Zm\'^  +  {i-QY\Vof^ 

771=1 

as  C?  — >  1  ,  (37) 

which,  as  expected,  is  the  same  case  found  in  equation  (23).  Secondly,  when  the  phase  variations 
become  completely  unpredictable  (i.e.,  p{(l>m)  =  1/(2^)  for  — tt  the  impact  is  that 

G  0.  In  that  case,  the  mean  of  the  beamforaier  output  power  spectrum,  shown  in  equation  (36), 
becomes 

M 

E{\V2\‘^}  (1  -  QY  E  <  ,  as  G  0  ,  (38) 

771=1 

which  is  the  power  output  of  an  unphased  (amplitude  summation)  array  with  random 
compensated  failures.  Finally,  when  there  are  no  element  failures  {Q  0)  but  there  are  still 
phase  perturbations,  the  mean  beamformer  output  power  spectrum  becomes 

M 

E{\V2\‘^}^{l-G^)f^widl{Q  +  G‘^\Vo\^,  as  Q ->  0 ,  (39) 

771=1 

which  is  the  mean  for  an  array  with  beamformer  output  spectrum 

M 

'^(^)  =  E  (40) 

771=1  » 

as  shown  by  Nuttall.^ 
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As  for  the  case  of  element  failure,  the  expected  value  for  the  case  of  random  phase  variations 
with  uncompensated  element  failures  follows  from  equation  (36),  and  is  derived  by  replacing 
{Q  —  1)^  (Q  —  1)  and  removing  the  cross-element  terms  from  the  element  summations  to  yield 

M 

m=l  (41) 

+  G^l-Qf  1^0^ 

Comparing  this  expression  with  that  in  equation  (36)  shows  the  effect  of  element  failure 
compensation  on  mean  beamformer  power  output  when  random  phase  variations  are  included. 

A  combination  of  different  physical  effects  may  cause  random  phase  errors.  In  that  case,  the 
net  phase  error  4>m  at  each  element  is  written  as  a  sum  of  the  phase  errors  (f>[^  from  the  various 
physical  effects.  Thus,  J  independent  physical  effects  causing  a  phase  error  produce 

<l>m  =  'Z4>'£,  (42) 

i=l 

which  results  in 

G  =  \E{jm}\, 

=  l■£^{exp(l]^</>^^)}|, 

=  (43) 

where  is  the  value  of  G  corresponding  to  the  source  of  the  phase  error.  This  expression 
follows  from  the  statistical  independence  of  the  various  sources  of  phase  error. 
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6.  RANDOM  AMPLITUDE  VARIATIONS 


The  final  stochastic  variation  under  consideration  is  random  amplitude.  The  derivation  begins 
by  building  upon  equation  (26)  through  the  addition  of  a  random  amplitude  variation  term  am- 
The  nominal  (no  amplitude  variation)  case  has  am  which  is  appropriate  because  no 
assumptions  have  yet  been  made  on  the  specific  statistical  nature  of  the  random  variables  in  this 
study  (i.e.,  they  need  not  have  zero  mean).  With  the  amplitude  variations  am^  the  beamformer 
output  spectrum  becomes 


M/2 

^^3(0)  =  E  /5m  {p^-mTiTn^m  “b  i  (44) 

771=1 

where  the  subscript  3  denotes  element  failures  plus  phase  and  amplitude  variations,  and  ^mi  1m 
are  defined  as  before.  In  the  analysis  that  follows,  the  only  assumption  made  for  the  random 
variable  am  is  that  it  is  statistically  independent  of  m.  Later  examined  is  the  special  case  where 
am  has  unit  mean,  which  represents  the  physically  important  example  of  unbiased  amplitude 
variations. 

As  in  the  previous  cases,  the  expected  value  of  the  beamformer  output  power  spectrum  is 
evaluated  as 


IVl  f 

£{inp}  =  E  (S{lAmin  -  |£{Am}P)  +  IBlV'alP. 

m=l 


(45) 


where 


Xm  —  Pmi^mlm^m  "b  (^mlm^m) 


and  a  dependence  on  u  is  assumed. 


The  first  term  in  expression  (45)  is  given  by 


=  E{Pl,}  ■  E{ai\'ym\'^\Zm\‘'  +  al\'rm\‘'\Zm? 

+  OimO^ilml^^m^m  “b  Iml'm^m^'m)}  j 

=  {l-Qf-{E{al}-i\Zm\^+\Zmf) 

+  E{amOim}  ■  |-^{77n}l  '  {^m^m  "b  > 

=  (I  -  Qf  ■  E{al}  ■  {\Zm\^  +  \Zmf) 

+  (1  -  Qf  ■  {E{am}f  ■  \E{jm}\^  ■  (z*^Zm  +  ZmZ*J  ; 


(46) 
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the  second  term  is  given  by 

=  \E{^rn}E{am}E{'ym}{Zm-¥  Zm)?  , 

=  \E{^m]?  ■  \E{am}?  •  ’  K  +  , 

=  (1  -  Q)^  ■  {E{am}f  ■  \E{jm}?  ■  \Zm  +  ;  (47) 

and  the  final  term  is  given  by 

l-E^Vs}!  =  Prni^mlfm^m  “h  Q^rnTm'^m) } I  j 

~  1  y  ^  ^{Pm\  ’  {E {(Xrn^m} •E'{^mTrn}'2^Th)  I  ? 

=  i(l-(3)^£{a„}B{7„}E(^«  +  ^)P> 

=  (l-(3)*-(£{a„})^-|£{7„}M^'oF-  (48) 

For  notational  convenience,  the  variable  Ak  is  defined  as  the  statistical  moment  of  am,  so  that 

Ak  =  E{a'^}  =  J  a’^p{am)dam ,  A:  >  0  ,  (49) 

which  is  independent  of  m.  Use  of  this  simplification  along  with  the  previously  defined  G  allows 
combination  of  the  above  terms  into  equation  (45)  and  (after  algebraic  simplification)  arrival  at 

E{\V,n  =  [A,{l-Qr-AlGHl-Qy]i:wldliCi) 

m=l 

W/2 

+  [yl?  G=  (1  -  0)2  -  Al  G2  (1  -  Q)‘  ]  E  (z„z4  +  2^2*)  (50) 

m=l 

+  AlG^l-QY  lUop. 

When  there  are  no  amplitude  variations,  cc^n  =  1,  so  that  Ai  =  ^2  =  1-  In  this  case,  the  mean 
beamformer  output  power  spectrum  of  equation  (50)  reduces  to  that  shown  in  equation  (36),  as 
expected.  Equation  (50)  represents  the  mean  beamformer  output  power  spectrum  for 
compensated  element  failures,  element  phase  errors,  and  element  amplitude  errors  combined.  It 
shows  the  important  interrelations  between  the  three  different  types  of  errors,  thus  justifying  a 
thorough  analysis  of  the  combined  effects. 

A  physically  important  class  of  amplitude  variations  is  the  unbiased  amplitude  variation, 
where  Ai  =  E{am}  =  1,  but  nothing  can  be  said  about  A2  =  E{a^'\.  In  this  case,  equation  (50) 
reduces  to 

M 

£{|VjP}  =  [A2{l-Q)^-G^l-QY]-£wi<PjQ 

m=l 
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iV±  I  ^ 

+  (1  -  Qf  -  (1  -  Qf]  x:  {ZmZl  + 


^G\1-QY\V,\\  (51) 

which  shows  the  separation  of  amplitude  and  phase  errors  (G  and  never  occur  in  the  same 
term). 

As  for  the  previous  cases,  the  results  are  compared  with  those  found  for  uncompensated 
failures.  As  before,  uncompensated  failures  are  represented  by  replacing  (1  -  QY  (1  -  Q)  and 
removing  the  cross-element  terms  from  the  element  summations,  which  yields 

M 

E{\V^,^ncomv?]  =  [A2(1-(5)-^?G2(1-Q)2]  E  <4(Ci) 

m=l  pzj 

+  AIGHI-Q?  l^oP 

for  the  general  case  with  uncompensated  failures  and 

m=l  pjj 

+  GH1-QY\Vo\^ 

for  the  case  of  unbiased  amplitude  errors  with  uncompensated  failures. 
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7.  SUMMARY  OF  VARIATIONS 


The  combined  effects  of  element  failures,  element  failure  compensation,  random  phase 
variations,  and  random  amplitude  variations  have  been  shown  in  the  preceding  sections.  In  all  of 
the  derived  expressions,  the  use  of  an  omnidirectional  sensor  directivity  response  and  no  error 
compensation  reduces  the  results  to  those  previously  reported.^’'*  The  detailed  derivations  have 
been  presented  to  illustrate  the  level  of  approximation  and  assumptions  involved  in  the 
development  of  the  mean  beampattems.  In  this  section,  the  principal  results  of  the  derivations  are 
gathered  for  easy  reference,  and  some  additional  comments  are  made  with  regard  to  their  use  m 
practice. 

The  two  primary  equations  used  for  array  error  analysis  are  equation  (52)  for  uncompensated 
failures  and  equation  (50)  for  compensated  failures,  both  of  which  are  repeated  here  for 
convenience.  For  the  general  case  of  uncompensated  failures,  the  mean  beamformer  output  power 

spectrum  is  given  by 

E{\Vuncomp\^}  =  [A2  (1  -  Q)  -  Af  (1  - 

+  AIGHI-QY 

and  for  the  general  case  of  compensated  failures,  the  mean  beamformer  output  power  spectrum  is 
given  by 

E{\Vcomp\^}  =  [A2  (1  —  (1  — 

+  [AlGHl-Q?-A\Cf‘{l-Q)*]Z  (55) 

*■  ^  771=1 

+  AIGH1-QY  IVoP, 

where  the  compensation  is  any  form  of  element  blanking  that  is  symmetric  to  failed  elements. 
The  second  summation  in  expression  (55)  has  been  simplified  from  the  previous  form  (see 
equation  (50))  using  the  simplified  relationship 

M/2  M 

771=1  771=M/2  +  1 

which  follows  from  the  definition  of  Zm-  The  parameters  Q  and  G  in  equations  (54)  and  (55) 
represent  the  failure  rate  and  the  amplitude  of  the  mean  error  due  to  phase  perturbations, 
respectively.  The  parameters  Ai  and  A2  represent  the  first  and  second  moments  of  amplitude 
errors,  respectively. 

In  these  expressions,  the  nominal  (no  error)  values  of  the  parameters  are  Q  =  0,  G  =  1, 

Ai  =  1,  and  A2  =  1-  Multiple  independent  phase  errors  with  independent  G’s  given  by  G^A  are 
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combined  to  have  G  =  Y[  Similarly,  multiple  independent  amplitude  errors  with 
independent  j4i’s  and  ^2’$  are  combined  to  have  >1:  =  0  and  ^2  =  11  ^2^^-  Nondirective 
sensor  elements  are  handled  by  modeling  the  sensor  directivity  as  omnidirectional;  hence, 
dm{Qi)  =  1  in  the  above  expressions  provides  equivalent  results  for  omnidirectional  sensors. 


The  first  term  in  equations  (54)  and  (55)  shows  that  there  is  a  sidelobe  floor  (minimal  sidelobe 
level)  due  to  errors  that  exists  regardless  of  the  level  of  sidelobe  control  for  the  nominal  array. 
Because  of  the  dependence  on  the  incident  signal  through  element  directivity  dm{Ci),  it  is  difficult 
to  quantify  this  effect  for  the  general  case.  However,  through  the  use  of  directive  array  gain, 
which  is  given  by 


AG  =  10  logic 


-  10  logic 


J 


an  approximate  sidelobe  floor  can  be  derived.  In  this  expression,  the  first  term  is  based  on  the 
element  response  in  the  steering  direction  Cs,  and  the  second  term  is  based  on  the  incident 
direction  Q.  Hence,  the  directive  array  gain  expression  AG  is  dependent  on  both  the  incident  and 
steered  directions. 

For  the  case  of  uncompensated  failures,  the  approximate  sidelobe  floor  is  given  (in  decibels) 
by 


min-(jS'Z/u7icojjip)-  ^  10  logic 


_ ^ _ 1 

AjGHl-Q) 


AG  . 


(56) 


As  expected,  improved  array  gain  lowers  the  level  of  this  sidelobe  floor,  even  in  the  presence  of 
errors  in  the  array.  For  the  case  of  compensated  failures,  the  level  of  approximation  applied  to  this 
sidelobe  analysis  is  similar  to  removing  the  cross  term  in  equation  (55),  in  which  case  the 
approximate  sidelobe  floor  is  given  by 


mm{SLcomp}  ~  10  logic 


A2 

AjG^l  -  QY 


AG  , 


(57) 


which  is  the  same  as  for  uncompensated  failures  with  the  working  element  rate  (1  —  Q)  replaced 
with  (1  —  QY-  Note  that  this  method  does  not  have  the  same  effect  as  doubling  the  failure  rate;  in 
general,  it  produces  a  very  small  difference  with  a  slightly  better  (lower)  sidelobe  level  for  the 
compensated  failures  as  compared  with  the  doubled  failure  rate.  This  behavior  is  the  result  of 
symmetry  that  is  maintained  because  of  the  particular  form  of  failure  compensation  employed. 

An  example  of  the  utility  of  these  sidelobe  floor  approximations  is  given  in  the  next  section. 
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8.  SAMPLE  ERRORS 


To  accurately  model  the  effects  of  various  errors  on  the  array,  the  impact  of  different 
combinations  of  Ai,  A2,  and  Q  needs  to  be  examined  with  regard  to  a  particular  array.  These 
are  the  fundamental  beamformer  error  parameters  for  beamformmg  an  array  of  directional 
sensors.  They  provide  the  basic  mathematical  mechanism  by  which  statistical  errors  enter  the 
expressions  for  the  statistical  beampattem.  In  this  set  of  parameters,  replaces  G  since  that  is 

the  only  form  in  which  G  enters  the  mean  beampattem  expressions. 

In  the  remaining  sections,  the  derivation  of  these  parameters  is  presented  for  common  types  of 
errors.  The  physical  mechanisms  that  affect  the  parameters  are  translated  into  effective  values  for 
beamformer  errors.  To  assess  the  effect  of  the  parameter  values  on  specific  beampattems,  the 
following  numerical  example  is  considered;  a  20-by-20  square  planar  array  of  elements  spaced  at 
Nyquist  sampling  (A/2  spacing  between  elements),  where  each  element  has  a  cosine  directivity 
pattern  in  the  direction  normal  to  the  array  surface.  The  array  is  shaded  in  both  directions  with 
Taylor  shading^  to  a  sidelobe  level  of  25  dB  with  n  =  4.  The  beampattem  for  this  array  under 
perfect  conditions  (no  errors  and  no  failures)  and  steered  to  the  array  maximum  response  axis  is 
given  in  figure  1  (only  a  slice  of  the  beampattem  through  one  direction  is  shown  due  to 
symmetry).  Figure  2  shows  the  effect  of  10%  and  20%  failures  on  the  beampattem  (Q  =  0.10  and 
Q  =  0.20,  respectively),  both  without  compensation  and  with  diagonal  blanking  compensation.  In 
figure  3,  the  effects  of  phase  errors  on  array  response  are  shown  for  values  of  G^  =  0.80, 0.90, 
and  0.95,  and  in  figure  4,  the  effects  of  unbiased  amplitude  errors  (Ai  =  1.00)  are  shown  for 
values  of  A2  =  1.05, 1.10, 1.15,  and  1.20.  Figure  5  illustrates  the  combined  effects  of  a  set  of 
phase  and  amplitude  errors  with  both  uncompensated  and  compensated  failures.  In  this  case,  the 
simulation  was  performed  with  G^  =  0.90,  Ai  =  1.00,  A2  =  1.10,  and  Q  =  0.10.  For  combined 
errors,  the  nulls  of  the  sidelobes  are  much  less  significant  than  for  any  of  the  component  errors. 
Thus,  the  combination  of  errors  both  increases  the  peak  sidelobe  level  (through  addition)  and 
smooths  out  the  oscillations  in  the  sidelobe  response.  This  smoothing  effect  is  more  significant 
with  large  errors  than  with  small  errors. 

The  utility  of  expressions  (56)  and  (57)  for  approximating  the  sidelobe  floor  can  be  shown  by 
considering  the  combined  error  case  in  figure  5.  For  uncompensated  failures,  the  sidelobe  floor 
expression  in  equation  (56)  was  used  to  obtain  the  approximate  floor  shown  in  figure  6,  which 
allows  an  assessment  of  the  improved  shading  in  the  array.  When  a  new  shading  scheme  is 
considered  to  lower  sidelobes,  two  effects  must  be  addressed.  The  first  is  that  the  floor  will 
decrease  according  to  the  increase  in  directive  array  gain  (AG).  However,  if  this  decrease  does 
not  bring  the  sidelobe  floor  below  the  level  of  desired  sidelobe  control,  the  shading  becomes  error 
dominated,  which  is  an  undesirable  condition  that  should  be  avoided.  If  conservative  error 
estimates  are  known  at  the  beamformer  design  stage,  the  shading  may  be  adjusted  to  manage  this 
sidelobe  control  versus  sidelobe  floor  tradeoff,  so  that  the  result  will  minimize  oscillations  in  the 
sidelobes  of  the  beampattem. 
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Incident  Angle  (deg) 


Figure  2.  Beampattern  for  the  Array  with  Failures 


Figure  3.  Beampattemfor  the  Array  with  Phase  Errors  of  Varying 


Figure  4.  Beampattemfor  the  Array  with  Amplitude  Errors  of  Varying  A2 


Incident  Angle  (deg) 


Figure  6.  Sidelobe  Floor  for  Combined  Errors  in  Uncompensated  Case 


8.2  ELEMENT  POSITION  ERRORS 

Element  position  errors  are  an  important  type  of  phase  error.  In  this  case,  the  individual 
element  phase  is  given  by 

<P^  =  {ki-ks)-Xm  (61) 

when  there  is  no  position  error.  It  is  assumed  that  the  beamformer  contains  no  knowledge  of  a 
positional  error  and  therefore  applies  phase  delays  according  to  the  assumed  (nominal)  position. 
Thus,  the  element  positional  errors  only  affect  the  ki  •  Xm  component  of  phase.  In  this  situation, 
the  phase  is  modified  to 

—  ki  •  {Xfn  +  Aa)  —  k^  •  Xm  —  {ki  ~  kg)  •  Xjn  +  ki  •  Ax  ,  (62) 

with  an  error  in  position  Ax. 

The  perturbation  in  this  case  consists  of  a  deterministic  vector  ki  (with  unspecified  direction 
and  known  magnitude)  dotted  with  a  random  positional  error  vector  (with  arbitrary  direction  and 
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random  magnitude).  The  error  in  position  can  therefore  be  stipulated  to  be  a  stochastic  magnitude 
with  an  arbitrary  direction,  resulting  in  two  independent  random  variables:  a  magnitude 
M  =  I  Ax  I  and  an  angle  6  that  is  the  angle  between  the  incident  wavevector  ki  and  the 
displacement  vector  Ax.  Thus,  the  dot  product  of  the  error  term  in  the  above  expression  can  be 
rewritten  as 


C"  =  -  rr  =  •  Ax  =  (cu/c)M  cos(^) .  (63) 

Because  the  error  in  displacement  vector  can  be  located  in  any  direction  (in  three-dimensional 
space)  with  equal  probability  and  the  direction  of  the  beamformer  steering  component  is  fixed  in  a 
statistical  sense  (i.e.,  it  is  deterministic),  the  combination  M  cos{9)  is  a  statistical  variable  that 
covers  all  angles  in  0  <  0  <  tt  equally  for  each  non-negative  value  of  M.  The  physical  meaning 
of  M  and  6  is  replaced  with  statistically  identical  variables  such  that  M  is  normally  (Gaussian) 
distributed  about  zero  and  9  is  uniformly  distributed  as  an  angle  over  the  upper  hemisphere  of  a 
unit  sphere  (it  is  arbitrarily  made  to  be  the  angle  relative  to  the  north  pole).  These  variables  still 
apply  in  the  same  way  to  but  the  use  of  a  normally  distributed  M  (allowing  values  less  than 
zero)  facilitates  the  solution.  An  alternative  approach  with  an  identical  solution  is  to  assume  M  to 
be  Rayleigh  distributed,  with  9  uniformly  distributed  as  an  angle  over  the  entire  surface  of  a  unit 
sphere. 

Because  the  angle  ^  is  a  random  variable  with  uniform  distribution  over  the  upper  hemisphere 
of  a  unit  sphere  (relative  to  the  north  pole),  it  has  a  cumulative  distribution  function 
c{9q)  =  Prob{9  <  9q)  given  by  (with  9  €  {0,7r/2}) 

\  r2-K  j-eo 

c(^o)  =  TT  I  /  cos(^)  d9  d(j)  =  sin(0o) ,  <  7r/2  .  (64) 

ZTT  Jo  Jo 

The  probability  density  function  of  the  random  variable  9  is  given  by 

d 

p{9)  =  -^c{9)  =  cos(^) ,  0  <  ^  <  7r/2  .  (65) 

The  random  variable  M  is  assumed  to  be  normally  (Gaussian)  distributed  with  zero  mean  and  a 
standard  deviation  given  by  cm-  Thus, 


is  the  probability  density  function  of  the  random  variable  M.  Note  that  this  variable  has  an  equal 
probability  of  being  positive  or  negative,  so  that  the  resulting  statistics  of  the  combination  of 
M  ■  q.os{9)  covers  a  statistically  uniform  sphere  with  random  magnitude.  The  distribution  of  the 
random  magnitude  is  given  by  the  right-half  (greater  than  zero)  side  of  the  normal  distribution  in 
equation  (66). 
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The  function  G  needed  for  the  mean  beamformer  output  power  spectrum  under  phase  errors  is 
evaluated  by  forming 


roo  /‘7r/2 


=  r  r\^’, 

J -CO  Jo 

1  roo  r7r/2 

J—oo  Jo 


p{e)p{M)dedM , 

/  iMu  cos{6) 


^  cos(9)dedM. 

2ah 


A  change  of  variables  from  M  to  nondimensional  M  =  (u)/c)M  and  from  0  to  nondimensional 
rj  =  cos(0)  is  used  to  reduce  the  double  integral  to 


/oo  rO  /  ^ 

J  exp  liMr] 


2a;2(7i 


dpdM , 


=  ^  ^  f 

OJy/^aM  -^0  \/l  —  '9^ 

Jo  y/l  —  T}^  \  2c^ 


MV 


ICJGm 


where  the  function  erf  (a;)  is  the  error  function  with  argument  x.  The  above  expression  uses  no 
analytical  assumptions  and  is  exact.  In  practice,  the  standard  deviation  cm  of  the  magnitude  of 
the  error  is  assumed  to  be  small  compared  with  a  wavelength.  Thus, 


which  implies 


Om  <  c//  , 


(7Mti;/(-\/2c)  <  e 


for  some  small  s.  The  assumption  of  small  cm  (compared  with  a  wavelength)  is  a  realistic  one, 
because  position  error  magnitudes  on  the  order  of  a  wavelength  would  render  a  phase-based 
acoustic  array  useless.  Since  the  argument  of  the  error  function  is  now  of  small  magnitude,  the 
small  argument  error  function  expansion  (derived  from  equation  7.1.5  in  Abramowitz  and 
Stegun®) 

9*  00  ^2n+l 


is  used  to  arrive  at  (after  some  algebra) 
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(70) 


where 


G  =  \E[^m\\  =  e-^  Z 


Sn 


^on\{2n+l)  ’ 


5  = 


2c2 


is  a  nondimensional  scaling  of  the  standard  deviation  of  position  error  magnitude.  In  practice,  the 
summation  in  this  expression  may  be  truncated  after  a  few  terms  (two  or  three  terms  is  usually 
enough)  for  small  position  errors. 


It  is  useful  to  compare  the  phase  error  due  to  positional  errors  with  the  highly  conservative 
approximation  to  phase  error  that  is  often  employed.  Restating  the  error  in  position  as  three 
independent  errors  in  orthogonal  directions  (say  x,  y,  and  z  for  convenience)  creates  a  phase  error 


of 


C"  =  «”  -  =  k!Ax‘  +  k^Ax’/  +  k!Ax‘,  (71) 

where  superscripts  are  used  to  represent  the  coordinate  directions.  The  random  positional  error  is 
now  modeled  with  each  component  (Ax^,  Ax^,  Ax^)  as  an  independent,  zero-mean  normally 
(Gaussian)  distributed  random  variable  with  standard  deviation  The  approach  followed  is  to 
bound  each  component  of  ki  by  its  magnitude  ko  =  u/c,  which  creates  a  conservative  bound  on 
the  phase  error  of 


IC"l<3«:o|Bl,  (72) 

where  D  is  a  random  variable  that  is  distributed  identically  to  the  components  of  Ax.  Using  the 
nominal  phase  error  result  of  the  previous  section,  G  is  now  given  by 

G  =  exp  =  exp(-95'),  (73) 

where  S  is  defined  as  before.  This  approach  is  identical  to  computing  a  bound  on  an  independent 
G  for  each  coordinate  direction  and  then  multiplying  these  bounds  together.  A  comparison  of  this 
very  conservative  bound  and  the  exact  term  (as  computed  in  equation  (70))  is  shown  in  figure  7. 

A  second  (and  more  accurate)  conservative  approach  is  developed  by  creating  the  phase  error 
in  the  form 


C"  =  kfAx=^  +  kfAxy  +  klAx^  =  fco  \D\  {k  ■  Ax) ,  (74) 

where  ki  and  Ax  are  the  unit  vectors  in  the  same  direction  as  ki  and  Ax,  respectively.  The 
dot-product  occurring  in  the  preceding  expression  is  now  conservatively  bounded  in  the  interval 
[—1,1],  such  that  the  phase  error  is  conservatively  estimated  by 
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Figure  7.  Phase  Error  (?  for  Highly  Conservative  Error  Approximation 


ICI  <  PI ,  (75) 

which  is  obviously  less  conservative  than  the  previous  bound.  This  estimate  leads  to 

G  =  exp  ^  2~^~^  =  exp (-5) ,  (76) 

with  the  same  parameter  S.  A  comparison  of  this  less  conservative  bound  and  the  exact  term  (as 
computed  in  equation  (70))  is  shown  in  figure  8. 

From  figures  7  and  8,  it  is  clear  that  the  addition  of  separate  bounds  on  the  different 
directional  components  of  error  causes  a  very  large  deviation  from  the  full  statistical  error.  This 
result  is  due  to  the  misrepresentation  of  ki  as  three  independent  vectors  of  equal  magnitude.  A 
better  conservative  bound  is  the  error  bound  in  equation  (75),  which  is  based  upon  a  bound  on  the 
dot-product  of  the  incident  wavevector  and  the  element  perturbation  vector.  Although  it  is  not 
recommended  to  use  this  bound  (since  a  more  exact  solution  in  the  form  of  equation  (70)  is 
already  available),  it  provides  some  practical  guidance  for  using  nonuniform  element 
perturbations.  If  the  perturbations  in  different  coordinate  directions  vary  in  magnitude  (such  as 
when  the  x-position  is  more  tightly  controlled  than  y  and  z),  then  a  reasonable  first  estimate  of  G 
can  be  obtained  by  taking 
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Figure  8.  Phase  Error  for  Less  Conservative  Error  Approximation 


Wm\  <  A:o-maxD,  (77) 

where  max  jD  is  a  random  variable  representing  the  deviation  in  the  dimension  with  the  maximum 
perturbation  (highest  standard  deviation).  This  approach  provides  better  estimates  than  the  one 
that  independently  adds  conservative  estimates  for  error  due  to  each  direction  (with  separate,  but 
identical,  on  each  term). 

The  preceding  method  uses  the  equivalent  phase  eixor  approach  outlined  earlier.  A  more 
accurate  solution  can  be  obtained  by  applying  the  perturbation  to  each  incident  angle  and 
independently  computing  the  resulting  beamformer  output  at  that  angle;  however,  that  level  of 
detailed  modeling  eliminates  any  insight  that  can  be  gained  from  the  first-order  effects  models. 


8.3  ELEMENT  ORIENTATION  ERRORS 


Element  orientation  errors  are  a  type  of  error  that  produces  closed-form  approximate 
statistical  results  in  specialized  cases.  For  the  analysis  here,  attention  is  restricted  to 
cosine-directive  elements  with  perfect  baffling.  In  this  case,  the  orientation  affects  the 
beampattem  through  the  element  directivity  term  dmiQ)  =  max{(^  •  rim,  0},  where  Q  is  the  unit 
vector  in  the  incident  wave  direction  and  rim  is  the  unit  vector  normal  to  the  element.  The  effect  of 
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the  element  orientation  error  is  a  resulting  random  perturbation  in  the  normal  vector  that  creates  a 
magnitude  change  in  the  element  directivity.  The  random  perturbation  can  be  represented  as 

dmiCi)  =  niax{5  •  (n^  +  An),  0}  =  max{(Ci  •  +  G  •  An),  0} ,  (78) 

where  An  is  the  perturbation  to  the  orientation  vector.  Since  the  incident  direction  Q  is  fixed  (in  a 
statistical  sense),  the  perturbation  term  Q  •  An  represents  a  random  magnitude  multiplied  by  the 
cosine  of  a  random  angle.  In  the  form  given,  the  perturbation  is  limited  to  the  constraint 
[iim  +  An|  =  1,  which  makes  the  connection  between  the  dot-products  and  the  angle  cosine 
meaningful. 

The  constraint  jn^  +  An]  =  1  has  the  following  geometrical  interpretation:  the  vector  n^  is 
given  as  a  unit  vector  (by  definition),  such  that  any  sum  of  this  vector  and  another  vector  that  is 
added  to  the  unit  vector  must  lie  on  the  unit  sphere  whose  origin  is  the  same  as  the  origin  of  rim- 
Some  geometrical  analyses  show  that  this  is  only  possible  when  the  added  vector  An  has  a 
magnitude  given  by 


|An|  =  2sin  ,  (79) 

where  4>m  is  the  angle  between  rim  and  n^  +  An.  The  variable  <j)m  now  has  a  physical 
interpretation  as  the  angular  error  in  orientation  for  the  element  of  the  array.  To  this  point,  the 
interpretation  of  |  An|  and  (j)m  is  independent  of  the  specific  element  directivity  of  the  array  being 
analyzed. 

Since  the  perturbation  Q  •  An  has  a  magnitude  given  by  |  An|,  the  only  remaining  component 
to  analyze  is  the  random  phase.  The  phase  of  Q  ■  An  is  analyzed  by  first  recognizing  that /or 
small  perturbations  /Isn  lies  nearly  on  the  plane  normal  to  the  vector  n^  (as  a  result  of  the  sum 
Hm  +  An  lying  on  the  same  unit  sphere  as  fim,  and  |  An|  being  small).  Defining  the  angle  ■0m  as 
the  angle  between  fim  and  Q  yields 


cos  7pm  =  fim 


The  angle  cp  is  defined  as  the  angle  between 


An  and  Q  such  that 


(80) 


coscp  = 


An  -  0 
An  I 


(81) 


A  final  angle  0  is  defined  as  the  rotation  of  the  perturbation  vector  An  in  the  plane  normal  to  the 
vector  fim,  so  that  0  is  the  random  variable  for  the  direction  of  the  orientation  error  for  the 
element.  Since  0  is  a  random  variable  with  uniform  distribution,  the  point  0  =  0  is  arbitrary  and 
can  be  set  to  be  the  location  where  the  projection  of  Ci  intersects  the  plane  normal  to  Um-  Some 
simple  geometry  shows  that  the  following  relationship  holds: 
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cos  V?  =  sin  ipm  cos  d . 


(82) 


Combining  equation  (79)  with  equations  (81)  and  (82)  yields 

Q  •  An  =  2  sin  sin  ipm  cos  6 


(83) 


as  an  expression  for  the  perturbation  under  small  orientation  errors.  Two  further  approximations 
simplify  this  formula.  First,  the  assumption  of  small  orientation  errors  implies  that  <  £ 
where  0  <  e  <§;  1,  so  that  2  sin(0,n/2)  «  Secondly,  directive  elements  whose  normal  is  far 
from  the  incident  wave  direction  contribute  little  to  the  beampattem,  which  implies  that  |V’m|  <  £ 
for  terms  that  contribute  significantly  to  the  beamformer  output.  Under  this  assumption, 
sin  ipm  ^  i’m  cos  iprn,  which  holds  for  the  elements  that  are  the  primary  contributors  to  the 
beampattem.  The  elements  that  are  the  most  significant  contributors  in  the  unperturbed  situation 
remain  the  most  significant  under  small  perturbation.  With  these  simplifications,  equation  (83) 
reduces  to 

Ct  •  An  «  4>mVm  cos  cos  9  , 

~  9')  •  Tlff^  .  (84) 

By  restricting  the  analysis  to  cosine  directive  elements  with  perfect  baffling,  equation  (78) 
becomes 

dm(Ct)  “  niaX'(Q:ni(Ci  ■  ^m))  O}  j 

OLm  ■  max{0  •  n^,  0} ,  (85) 

where  the  approximation  holds  for  significant  contributors  to  the  beampattem  (elements  with 
nontrivial  directivity)  and 


Q!jn  =  1  +  COS  9 


(86) 


is  used  to  represent  the  element  magnitude  perturbation  due  to  orientation  errors  for  cosine 
directive  elements. 

To  use  the  statistical  analysis  developed  in  earlier  sections,  the  random  variables 
(comprised  of  9,  and  ’tpm)  must  be  independent  and  identically  distributed.  The  variable  (f)m 
represents  orientation  error  magnitude  (in  radians),  or  tilt,  and  the  variable  9  is  the  direction  of 
this  tilt,  so  that  these  variables  are  assumed  to  be  identically  distributed.  However,  the  variable 
V’m  represents  the  angle  between  the  incident  vector  and  the  vector  normal  to  the  unperturbed 
element,  which  is  not  random  and  necessarily  varies  among  the  elements.  Since  this  variable  must 
be  identical  for  all  elements,  an  approximation  must  be  used.  Two  approximations  have  been 
tested:  one  sets  'tpm  to  be  the  average  of  all  for  each  incident  angle  and  the  other  sets  ipm  to 
the  maximum  'ipjn  for  each  incident  angle.  Because  the  former  is  easier  to  implement  and  the 
difference  in  the  two  was  seen  to  be  minor. 
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V’m  ~  S  arccos(Ci  •  nj) ,  •  Vm  ,  (87) 

is  a  reasonable  approximation  for  arrays  with  small  curvature.  Under  this  restriction,  the  problem 
is  appropriate  for  the  analysis  developed  in  earlier  sections. 

The  physical  interpretation  of  the  two  random  variables  implies  that  6  is  uniformly  distributed 
over  {-• TT,  tt}  and  4>m  follows  a  distribution  similar  in  shape  to  the  right  side  (positive  component 
only)  of  a  normal  (Gaussian)  distribution  with  zero  mean.  To  facilitate  the  analytical 
simplifications,  the  random  variable  9  is  restricted  to  {0,  tt}  and  (j)m  is  allowed  to  follow  a  true 
Gaussian  distribution  (permitting  both  positive  and  negative  values)  with  zero  mean  and  a  stated 
standard  deviation  of  ag.  In  this  way,  the  variable  (j)m  is  allowed  to  have  values  ranging  from 
{-00,  oo}  even  though  the  values  outside  of  {-tt,  tt}  are  nonphysical.  But  for  small 
perturbations,  these  values  are  on  the  extreme  tails  of  the  distribution  and  therefore  are  negligible. 
Thus,  (To  is  the  standard  deviation  of  element  orientation  (tilt)  with  respect  to  normal,  given  in 
radians,  and  ^ 

The  necessary  statistics  on  am  are  given  by 
Ai  =  E[am] , 

/4-00  rTT 

/  {I +(t>m'(pm  cos  6)  p{(j)m)p{0)ded<{)m, 

-00  Jo 

=  - K=  r°°  [  {I  +  COS  6)  exp  (  )  dd  dcpm , 


which  shows  the  perturbations  to  be  unbiased  (i.e.,  equation  (51)  holds),  and 
M  =  E\o?^ , 

/-hoc  pTT 

/  {\ +<j)m'4’m  cos  efp{4)m)p{d)ded(i>m, 

-00  •/O 


1 

r+oo  rir 

'  /  (1 +  (^„i'0mCOS^)^exp 

“00  Jo 

TraoV^  J 

_  1  / 

(2  +  (f>m^l)exp  ' 

I  j  T  TTL  0 

“  2 

*  d(j>Tn  5 


for  which  equation  (87)  is  used  to  obtain  the  value  of  'tpm  for  each  incident  angle. 
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8.4  ELEMENT  RESPONSE  TOLERANCES 


Element  response  tolerances  are  a  simple  form  of  error  that  can  be  modeled  by  placing 
statistical  variations  on  the  magnitude  of  individual  sensors.  If  these  variations  are  assumed  to  be 
normally  distributed  about  the  mean,  a  simple  normal  (Gaussian)  distribution  is 


p{am)  =  — ^  exp 

(7aV27r  V  2al 


vi^here  0;,^  is  the  magnitude  of  the  sensor.  In  this  statistical  model,  it  is  assumed  that  the 
nominal  (mean)  sensor  response  is  am  =  1,  and  the  standard  deviation  from  this  mean  is  given  by 
(7a.  If  element  response  tolerances  are  given  as  a  standard  deviation  in  received  power  (in 
decibels),  then  a  transformation  between  Om  and  Pm  (the  tolerance  in  decibels)  shows 


dm  = 


where  Pm  is  now  a  random  variable  with  a  mean  of  zero  and  a  standard  deviation  of  a^,  which  is 
the  stated  standard  deviation  in  power  level.  For  this  case,  the  errors  may  have  a  bias  in 
magnitude;  thus,  the  parameter  Ai  must  be  computed  as 


Ai  =  E{am}  =  I  dp , 


which  is  best  evaluated  numerically.  In  a  similar  manner,  the  parameter  A2  is  given  by 

which  is  also  best  evaluated  numerically. 

The  use  of  a  normally  distributed  magnitude  variation  as  in  equation  (91)  can  be  easily 
exarmned  analytically.  In  this  case,  the  required  statistical  moments  Ai  and  A2,  as  found  in 
equation  (50),  are  given  by 

Ai  =  E{am}  =  1 


E{al} , 

E{{am-pa?}  +  pl, 


^a  +  1: 
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where  Ha  =  E{am}  =  =  1  is  the  mean  value  and  the  third  step  follows  from  the  definition  of 

variance.  These  expressions  can  then  be  used  in  equation  (50)  to  obtain  the  effects  of  element 
response  tolerances  on  mean  beamformer  spectral  power  output. 
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9.  SUMMARY 


The  error  analysis  of  a  two-dimensional  passive  array  (conformal  or  planar)  with  directive 
element  responses  has  been  developed.  The  analysis  was  based  on  a  set  of  fundamental 
beamformer  error  parameters  that  address  effective  element-level  phase  errors,  amplitude  errors, 
and  statistical  failures  (reliability).  The  combination  of  these  various  parameters  and  their  effect 
on  mean  beampattem  power  output  was  given.  Also  presented  were  numerical  examples  of  the 
errors,  as  well  as  the  analytical  derivation  of  various  physical  error  mechanisms  in  terms  of 
amplitude  and  phase  error.  Because  the  mathematical  derivations  described  here  pertain  to  a  wide 
variety  of  passive  arrays  (and  include  the  effects  of  array  curvature  and  element  directivity  with 
the  errors),  these  expressions  will  have  broad  application  for  assessing  the  impact  of 
element-level  errors  on  array  beamformer  performance. 
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